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Abstract. Inspired by Quantum Mechanics, we reformulate Hilbert's tenth problem in the 
domain of integer arithmetics into problems involving either a set of infinitely-coupled non-linear 
differential equations or a class of linear Schrodinger equations with some appropriate time- 
dependent Hamiltonians. We then raise the questions whether these two classes of differential 
equations are computable or not in some computation models of computable analysis. These are 
non-trivial and important questions given that: (i) not all computation models of computable 
analysis are equivalent, unlike the case with classical recursion theory; (ii) and not all models 
necessarily and inevitably reduce computability of real functions to discrete computations on 
Turing machines. However unlikely the positive answers to our computability questions, their 
existence should deserve special attention and be satisfactorily settled since such positive answers 
may also have interesting logical consequence back in the classical recursion theory for the 
Church- Turing thesis. 

Monday, July 4, 2005 



Introduction 

Hilbert's tenth problem [21 E] is concerned with the availability of a universal procedure to 
determine whether an arbitrarily given Diophantine equation 

(1) D{xu---,xk) - 

has any positive integer solution or not. After more than 70 years since its inception, it has finally 
been shown that the problem is recursively noncomputable since no such universal procedure which 
is also classically recursive can exist. The problem is equivalent to the Turing halting problem and 
intimately links to the concept of effective computability as defined by the Church- Turing thesis 
in classical recursion theory. 

Nevertheless, we have established elsewhere, through an inspiration provided by quantum me- 
chanics 01 El 13 El El I some surprising connections of the above problem in number theory to 
problems over the continuous variables. In particular, we have been able to reformulate Hilbert's 
tenth problem in terms of a set of infinitely coupled non-linear differential equations for any given 
Diophantine equation |S1. Also, through the framework of quantum adiabatic computation [31, 
we have also associated Hilbert's tenth problem with a class of linear Schrodinger equations with 
appropriate time-dependent Hamiltonians. It has been proposed then that a physical implementa- 
tion of quantum mechanical processes for these Schrodinger equations could provide the physical 
means to solve Hilbert's tenth problem Here, in this paper we will only concern ourselves 
with the differential equations as mathematical objects, however, and will not appeal to any real 
physical processes. The mathematical objects so derived are extremely valuable as they provide 
us a direct link between the well-established classical recursion theory and the infantile subject of 
computable analysis, through a known noncomputable problem in the former theory. 



2 



TIEN DUNG KIEU 



We especially want to raise in this paper the computability questions for these differential 
equations in the domain of computable analysis, outside and encompassing the classical recursion 
theory where Hilbert's tenth problem was originally formulated. 

From the recursive noncomputability of Hilbert's tenth problem, one might conclude that these 
differential equations should also be noncomputable in the wider framework. Such hasty conclu- 
sion, however, is not warranted because of several reasons. Firstly, however likely the case one 
might expect, it would need to be established rigorously as a mathematical truth -because the un- 
solvability of Hilbert's tenth problem is only established in the framework of Turing computability, 
not necessarily in mathematics in general. 

Secondly, there are many computation models in computable analysis |2 IIOL I13L I14L ITB) but 
they are not all equivalent. And it is known that computability in one model may not be the 
same in some other model, see jl5j . for example, for a brief discussion comparison of the various 
models. This situation is in stark contrast to the classical recursion theory of functions from 
N to N. There, many different formulations have been given (notably by Kleene, Turing, Post, 
Herbrand/Godel, Markov) but in the end these all lead to the same notion of computability with 
the same class of computable functions. Such an equivalence has led to the postulation and support 
of the Church- Turing thesis in that theory. 

The computability definitions for a single real number in most, if not all, different computa- 
tion models of computable analysis are equivalent. However, for sequences of reals the definitions 
diverge and are not all equivalent. (The discussions on sequences of real numbers are necessary 
because of the necessity of topological notions in analysis.) This divergence results in the depen- 
dence of the notion of computability on the different computation models,^ or even on different 
choices within a single model. ^ 

In view of such an inequivalence of computability in different approaches, it is not at all a forgone 
conclusion that the noncomputablity of the differential equations mentioned above is trivially the 
only, inevitable possibility. 

In the next Section we brieffy present the observation inspired by quantum mechanics that leads 
to the connections between Hilbert's tenth problem with unbounded self-adjoint operators acting 
on some infinite-dimensional Hilbert space. From this we then reformulate Hilbert's tenth problem 
in terms of a set of infinitely coupled non- linear differential equations (eqs. (|15l below). We 
also propose a so-called continuation procedure to approximate some relevant part of its solution; 
the computability of the continuation procedure is left as an unanswered question for the time 
being. 

We then present a class of linear Schrodinger equations (eq. H18|l below) with a special class of 
time-dependent Hamitonians that are intimately connected to Hilbert's tenth problem. We call 
this the dynamical approach^ because of its use of the Schrodinger equations, to distinguish it from 
the kinematic approach above from which we derive the set of non-linear differential equations. 
We conclude the paper with some remarks and a discussion on the possible implications, including 
those on the Church- Turing thesis, if the differential equations are indeed computable in some 
computation model of computable analysis. 



^For example, the so-called real-RAM approach [nil4| introduces computable real functions directly and borrows 
only the concept of control structure from Turing computation, without any further referencing to the latter. The 
computability notion in this model, as a result, is different from that of other real computation models that are 
based and built from the Turing computation. 

■^One famous example is that the choice of different norms in a Banach space can lead to opposing conclusions 
about the computability of solutions of the same wave equation in three spatial dimensions with computable initial 
functions |13|. 
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HILBERT'S tenth PROBLEM AND UNBOUNDED OPERATORS IN HiLBERT SPACES 

Given a Diophantine equation with K unknowns x's as in eq. we can make a connection, 
following 01 1^1 El ISl IHl , with the following self-adjoint operator acting on some appropriate Fock 
space (a special type of Hilbert space) 

(2) Hp ^ {p{a\ai,--- ,a^j^aK)f , 
where 

(3) [aj,a\]^5jk, [aA;,aj]=0, 

which are usually termed the creation and annihilation operators, and most commonly seen in 
text-book treatment of the quantum simple harmonic oscillators. The Fock space is built out of 
the "vacuum" by repeating applications of the creation operators a]. 

The operator has non-negative and discrete eigenvalues (Z3(ni,--- ,rtK))^, with natural 
numbers ni, • • • ^riK- There is an eigenstate \Eg) corresponding to the smallest eigenvalue Eg.^ 
If the self-adjoint operator is considered as a Hamiltonian for some dynamical process then these 
are respectively the ground state and its energy. 

It is then clear that the Diophantine equation ^}) has at least one integer solution if and only 

if Eg = (^D (^n''i\ ■ ■ ■ ,n^^^^^ — 0, for some K -tuple of natural numbers (^n^\ ■ ■ ■ ,n^x^ ■ 

To sort out this Eg among the infinitely many eigenvalues is almost an impossible task. The 
strategy we will employ, as inspired by quantum adiabatic processes, is to tag the state \Eg) by 
some other known state \Ei) which is the ground state of some other self-adjoint operator Hj, 
which can be smoothly deformed to Hp through some continuous parameter s G [0, 1]. To that 
end, we consider the interpolating operator 

= Hi + f{s){Hp-Hi), 

(4) = Hi + f{s)W, 

which has an eigenproblem at each instant s, 

(5) [^{s)~Eq{s)]\E,{s)) ={), Q = 0,1,--- 

with the subscript ordering according to the sizes of the eigenvalues, and /(s) some continuous 
and monotonically increasing function in [0, 1] 

(6) /(0) = 0; /(1) = 1. 

Clearly, Eo{0) = Ej and Eo{l) = Eg. 
A suitable Hj is, where a's e C, 

K 

(7) = X! '^'("i " - "i): 

i=l 

which we will employ from now on. Here, Ej — and \Ej) = \ai ■ ■ ■ ax) is the Cartesian product 
of the coherent states 

(8) |«,) = e-^Y. 

ni=0 



■^Assuming that wc could always, if we need to, eliminate any degeneracy of the eigenvalue by, for example, 
modifying the Diophantine equation or by adding some small perturbation terms in Hp. 
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where \ni) are the eigenstates of ajai with eigenvalues Ui. The Xi can be chosen to be rational or 
even irrational numbers such that the first order equation 

K 

(9) ^A,K = 0, 

has no integer solutions in pi. This condition is to ensure that all the eigenvalues of Hi are non- 
degenerate, since the eigenvalues of Hj, which are of the form X^iLi -^i^-ij then easily seen to 
be unique for different if-tuples of natural numbers n's. 

The spectral flow - The "kinematic" approach 

We now derive the differential equations for the tagging connection between the instantaneous 
eigenvalues and eigenvectors at different instant s in |SJ- 

Note firstly that, from the normalisation condition {Eq\Eq) = 1, we can write 

for some real (j)q . This can be absorbed away with the redefinition 

(10) e^<^'(^)|£;,(.)) ^ 

upon which 

(11) {E,\ds\Eq) = 0. 
Differentiating ^ with respect to s yields 

(12) [f'{s)W-d,Eq]\Eq) + [f)-Eq]ds\Eq) = 0. 
We next insert the resolution of unity at each instant s, 



1 = ^ |i?™(s))(i?™(s)|. 



just after Sj in ifT^ to get, by virtue of |[TT|l . 

oc 

(f3) Egds\Eg) = [f'{s)W--dsE,]\E,) + Y,E„^{E,n\ds\E,)\E^). 

The inner product of the last equation with \Ei) gives 

(14) {Eg - Ei){Ei\d,\E,) = f{s){Ei\W\Eg)-d,E,5gi. 

Thus, for g ^ Z this gives the components of ds\Eq) in \Ei), provided Eq{s) ^ Ei{s) at any 
s e (0, 1), a condition whose proof has been outlined in 6 . Consequently, together with 

(15) dAE,) . /'(.)f:»^ii^,). 

Also, putting q = I in 114|) we have 

(16) d.Eqis) = f{.s){Eqis)\W\Eqis)). 

Equations (|15|l and 1)16(1 form the set of infinitely coupled differential equations providing the 
tagging linkage we have been looking for. 

In this reformulation, the Diophantinc equation ^ has an integer solution if and only if 

(17) lim£;o(s) = 0, 

s— »1 
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from the constructively known eigenvalues and eigenstates of Hj as the initial conditions. The 
limiting process might be necessary since Hp, i.e. may have a degenerate ground state 

eigenvalue in general. 



The non-linear differential equations of the last Section are infinitely coupled and may not 
be solved explicitly or computably in general. However, we are only interested in the ground 
state eigenvalue Eq{\) being zero or not. And since the influence on the ground state by states 
having larger and larger indices diminishes more and more thanks to the denominators in (|15|l . 
this information might be derived in some approximation scheme in which the number of states 
involved is truncated to a finite number. The size of the truncation cannot be universal and must 
of course depend on the particular Diophantine equation under consideration. 

In the below we speculate on an analytic approximation under the name of continuation pro- 
cedure, and we make no claim about its computability here but leave it as a challenge for the 
future. 

• Starting from the initial condition comprising of the constructively known eigenvalues and 
eigenvectors of 7J/ at s = 0, the differential equations ifT^ and [TE|l give us the series 



• If the remainders and 0i above are computable, we would be able to evaluate the 
radii of convergence in s, which contain s = ei, and also the truncation size Ni which 
determines the accuracy of the expansions.. 

• We then proceed to evaluate new series approximations similar to the ones above but this 
time centred at s = ei . The new series have new radii of convergence and a new truncation 
of N2 eigenvectors previously approximated at s = ei. After this step we have then covered 
a finite domain in s away from zero, with some computable degree of accuracy. 

• We keep reiterating this procedure, if possible, to obtain new remainders and radii of 
convergence and thus extend the covered domain in s, until we could evaluate the limit 
lims^o Eq{s). 

This is reminiscent of the procedure of analytic continuation of functions in complex analysis. 

Note that at each step of the above procedure we only require some finite truncation Ni for a 
given accuracy of the series approximations. (That accuracy would also "recursively" determine 
the truncations Ni at all previous steps, j < i.) Note also that, in general, the condition a; = 
for a computable real number x may not be effectively decidable in some computation model. 
But here we have the imposed condition that the eigenvalues at s = 1 must be integer-valued. 
This additional condition might help making the equality condition effectively decidable at s = 1 . 
That is, we would only need to approximate Eq(1) by the procedure above up to some accuracy, 
say 0.3, which sufficiently enables us to distinguish different integers, and we would not require 
infinite precision. The imposed condition of integer- valued eigenvalues for Hp, by construction, 
would provide us the built-in infinite precision at no extra cost! 

HILBERT'S TENTH AND THE SCHRODINGER EQUATION - THE "DYNAMICAL" APPROACH 

The decision result for Hilbert's tenth problem can also be encoded in yet another class of 
linear differential equations, apart from the class of nonlinear equations H15II16|I above. The linear 



Kinematic" continuation procedure? 



expansions 





E,{Q) + €^fmE,{{))\W\Eq{Q)) + Qi. 
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equation is just the Schrodinger equation which captures the dynamics of a proposed quantum 
adiabatic algorithm [2 El 13 El ■ 

Let li^it)) be the quantum state at time t (of some quantum system), its time evolution is given 
by the Schrodinger equation, for < t < r, 

(18) dtim) = ~iHt/r)\m), 

IV'(O)) = \ai---aK), 

where we have chosen the initial state at time t = to be the non-degenerate ground state of Hj. 

Once again we are only interested in the ground state of i5(l) — Hp. The quantum adiabatic 
theorem ^2] asserts that as r — + cx3 (that is, when the Hamiltonian Sj{t/T) in eq. (|18|l varies 
sufficiently slowly in the time t) the state of the above system at i = r would be in the desired 
ground state, IV'C''")) ~ l^g): to any arbitrary degree of precision! For a given precision, various 
versions of the theorem dictate different conditions on r in terms of some intrinsic properties of 
the eigenvalues and eigenfunctions of ^{s). Those conditions thus are highly dependent on the 
individual Diophantine equation, and hence are not suitable to the spirit of a universal procedure 
required by Hilbert's tenth problem. 

We have proposed a different and universal criterion for the identification of the ground state 

\Eg) from \iP{t)): the Fock state n^i\ - • ■ ground state \Eg) if it has an occupation 

probability greater than one-half. That is. 



(19) |(^(r)K°),...,n(;!))|'>l/2 |<\ • • • , n^^^) = ji?,). 

It should be emphasised here that such a criterion should be taken at the present time as a 
postulate as it has only been proved in some limited settings [51 . With this criterion, we only need 
to repeatedly solve the Schrodinger H18() for larger and larger r each time until the probability 
condition is satisfied so that the ground state can be identified. Such a time r can be shown to 
exist and be finite by the quantum adiabatic theorem. More details of these can be found in [HlEl- 

Note also that while the criterion H19|l may not the only one suitable for a physical imple- 
mentation of the Schrodinger equation, it is the only one that we could yet find suitable for the 
mathematical discussion of this paper. 

Unlike the case of nonlinear differential equations of a previous Section, here we could try to 
make use of a powerful computability result in a computation model which is known as the First 
Main Theorem by Pour-El & Richards ^Sl- Essentially, the theorem asserts that a bounded linear 
operator from a Banach space to a Banach space which maps a computable sequence of spanning 
vectors into another computable sequence will also map any computable element into another 
computable element. For the case at hand, our Schrodinger equation defines a linear operator, 

(20) U{t) ^ Texp|-ir^ Sj{s)ds 

where 1 is the time-ordering symbol, which maps the initial state to the final state in the same 
separable Hilbert space. Now, our initial state \ai ■ ■ ■ ax) is computable by construction. On the 
other hand, the linear operator H2Ufl coming out of the Schrodinger equation must be unitary and 
thus be bounded. Hence, the conditions of the theorem remained to be checked for fulfillment are: 
(i) which mathematical norm should be chosen to enable the identification criterion p9|l . or some 
other equivalent criterion, for the ground state at t = r; and (ii) whether the image of a particular 
computable basis is computable or not with this suitably chosen norm. 
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Concluding remarks 

Inspired by quantum mechanics, we have reformulated the question of solution existence of a 
Diophantine equation into the question of certain properties conceived in an infinitely coupled set 
of nonlinear differential equations. In words, we encode the answer of the former question into the 
smallest eigenvalue and corresponding eigenvector of a self-adjoint operator whose integer-valued 
eigenvalues are bounded from below. To find these eigen-properties we next deform the operator 
continuously to another self-adjoint operator whose spectrum is known. Once the deformation is 
also expressible in the form of a set of nonlinearly coupled differential equations, we could now 
start from the constructive knowns as a handle to study the desired unknowns. 

In addition, we also explicitly present a class of linear Schrodinger equations whose solutions 
at some time t from appropriate initial conditions contain the decision results for any given 
Diophantine equation. 

These reformulations map a noncomputable problem in the domain of integer arithmetics into 
the wider framework of computable analysis. We have given the names, as explained in the 
Introduction, "dynamical" and "kinematic" respectively for the resulting linear and non-linear 
equations. In particular, our set of nonlinear equations would be an important topic for the 
largely untouched subject of nonlinear computable analysis. 

For the various reasons given in the Introduction Section, the questions of computability for 
these differential equations are non-trivial and important. Towards some answers for these ques- 
tions, we have advocated and speculated on an approach based on the First Main Theorem by 
Pour-El & Richards for the Schrodinger equations, and some other approximation procedure 
for the set of nonlinear differential equations. However, for now, these computability questions 
will have to be left as open problems and challenges. 

If in the (admittedly unlikely) event that there exists computation model (with suitably cho- 
sen norm) in which the differential equations above are computable and that this model could be 
restricted and applied to functions from N to N then Hilbert's tenth problem would seem to be 
solvable in integer arithmetics (as opposed to be solvable in, for example, some physical quan- 
tum adiabatic computation)! And this would entail a logical breakdown of the Church- Turing 
thesis. Because the restriction of such a computable analysis model to the domain of integers 
would provide a new notion of effective computability different from that of the class of Turing 
computation. 

To be sure that this is not the case, further investigations are urgently needed to find out (i) 
whether the above equations are noncomputable in a// computation models in computable analysis; 
or (ii) whether the model that admits such computability, if exists, cannot be restricted to integer 
arithmetics. 
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